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具 Monod-Haldane 功 能 反应 的 四 种 群 食 钮 -捕食 
脉冲 系统 的 动力 学 分 析 * 


RLE, RMA, F R? 
(1- 桂林 电子 科技 大 学 数学 与 计算 科学 学 院 ， 桂 林 541004; 2- 合肥 师范 学 院 数学 系 ， 合 肥 230061) 


摘 ”要 : 基于 害虫 综合 管理 策略 ， 本 文 利 用 脉冲 比较 定理 、Floquent 理论 及 微小 扰动 法 ， 研 究 了 一 类 具 
有 Monod-Haldane 功能 反应 、 脉 冲 比例 收获 和 脉冲 常数 投放 的 四 种 群 食 饵 -捕食 系统 的 动力 学 性 
质 ， 证 明了 系统 两 食 饵 灭绝 和 持续 生存 的 充分 条 件 ， 而 且 给 出 了 一 食 饵 种 群 灭 绝 其 余 三 种 群 持续 
生存 的 两 个 充分 条 件 。 数 值 模拟 表明 ， 随 着 投放 量 的 增加 ， 系 统 出 现 拟 周期 分 支 、 倍 周期 分 支 、 
混沌 、 半 周期 分 支 等 复杂 的 动力 学 性 质 。 
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1 ”系统 的 提出 


人 类 与 有 害 的 生命 体 的 斗争 进行 了 几 千 年 ， 有 害 昆 虫 和 某 些 有 害 带菌 者 对 重要 的 植物 、 动 
物 和 人 类 造成 危害 ， 如 何 使 这 种 危害 减少 到 最 低 一 直 是 昆虫 学 家 与 社会 关心 的 问题 。 近 几 年 ， 
脉冲 微分 方程 四 被 引入 到 种 群 动力 学 中 ， 并 且 得 到 了 广泛 的 应 用 2 了 。 对 生物 种 群 的 控制 主要 
有 三 种 最 常见 的 方法 且 各 有 其 效 ， 即 化 学 防治 、 生 物 防治 图 MRS AP, HpaI 
略 是 使 经 济 损失 减 到 最 小 的 一 种 非常 有 效 的 方法 。 

食 饵 -捕食 系统 是 一 种 非常 重要 的 生物 系统 ， 已 有 许多 学 者 对 此 做 了 大 量 的 工作 03 区 ， 近 
几 年 的 工作 主要 集中 在 对 具有 功能 反应 的 食 饵 -捕食 系统 的 研究 上 。 对 具有 Monod-Haldne 功能 
反应 0el 和 脉冲 效应 食 饵 -捕食 生物 模型 的 研究 亦 有 许多 结果 08748]， 但 是 采取 综合 防治 策略 ， 
且 脉 冲 投放 发 生 在 不 同 国定 时 刻 的 结果 很 少 。 文 献 [19] 研究 了 四 种 群 食物 链 方 程 的 全 局 渐 近 稳 
定性 ， 向 中 义 、 宋 新 宇 等 学 者 对 一 类 具有 脉冲 效应 且 脉 冲 发 生 在 一 个 固定 时 刻 的 四 种 群 线性 食 
饵 -捕食 系统 作 了 分 析 P0， 随 后 又 以 相似 的 方法 讨论 了 脉冲 发 生 在 两 个 固定 时 刻 同 一 个 线性 模 
型 Ei。 仅 依 靠 天 敌 ， 有 了 时 不 足以 消灭 所 有 害虫 ， 喷 酒 杀 虫 剂 的 过 程 中 往往 也 会 伤害 天 敌 。 如 
果 根 据 食 饵 与 捕食 者 种 群 的 不 同 特性 ， 采 取 在 不 同 的 固定 脉冲 时 刻 ， 分 别 对 所 有 种 群 进行 收 
获 (消灭 ) 和 引进 天 敌 的 办 法 ， 就 既 能 起 到 对 天 敌 的 保护 ， 又 能 有 效 的 杀 灭 害虫 ， 这 样 更 能 有 效 
地 实现 对 系统 的 控制 ， 而 且 比 脉冲 发 生 在 一 个 固 时 刻 模 型 更 能 准确 地 描述 实际 问题 。 当 然 ， 我 
们 的 主要 目的 在 于 控制 害虫 生长 而 又 保护 天 敌 的 生存 ， 实 现 生 态 系 统 的 平衡 。 

本 文 研究 下 面 一 类 复杂 的 具有 Monod-Haldne 功能 反应 、 脉 冲 比例 收获 (消灭 ) 和 脉冲 常数 
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投放 发 生 在 不 同 固定 时 刻 的 四 种 群 食 饵 -捕食 者 系统 ， 并 讨论 其 全 面 的 动力 学 性 质 。 
ai(t) = zı (t)(dı — z1 (t) — azxe(t)) 
arzı (t)yı (t) a271(t)y2(t) 
1+0121(t) + cr?(t) 14 O2r1(t) + c2z2(t)’ 
t2(t) = x2(t)(d2 — Bxi(t) — x2(t)) 
bi z2(é)y1 (t) baxa (t)y2(t) PARETE DI tni; 
1+ pize(t) + pixk(t) 1+ pexe(t) + poad(t)’ 


; = yaxı (t) yb1x2(t) 

i (t) ob 人 Een eae eT STEP aap) 
nazz (t) nb2x2(t) ) 

1+ O2m1(t) + c22? (t) — 1+ poxe(t) + porz (t) 























a(t) = y2(t)( — da + 

A _ (1) 
xi(t) = —qi21(t), 

Ax2(t) = —qox2(t), 

Ayi(t) = —q3y1 (t), 


Ayo(t) = —qaya(t), 
Azı (t) = 0, 


t= (n+!l~ 1)7, 


Ax2(t) = 0, 

AYi(t) = pi, 

人 Ay2(t) = p2, 
其 中 z(t#) (i = 1,2， 以 下 表示 相同 ) 表示 两 食 饵 种 群 在 1 时 刻 的 密度 ，y(t) 表示 两 捕食 者 种 群 
在 t 时 刻 的 密度 ， 

p(z1) = aya; /(1 + O71 + c12?) 

称 为 Monod-Haldane 功 能 反应 函数 ,了 为 脉冲 周期 ，d; > 0 为 两 食 钮 种 群 z(t) 的 内 裹 增长 
E, dip 为 捕食 者 yi(t) HECE, a, 8 分 别 表示 两 食 饵 z(t) 种 群 相互 作用 对 其 相应 增长 率 的 
影响 ，a; > 0, bi > 0 是 食 饵 种 群 因 捕获 而 减少 的 比例 ，0 < oj < 1(7 = 1,2,3,4) 为 (n 填 1 一 
1)T 了 时 刻 对 四 种 群 的 捕获 或 消灭 比例 ，p; > 0 为 脉冲 大 小 或 迁 入 量 ， 表 示 由 于 人 工 繁殖 或 释放 
捕食 者 在 nT 时 的 密度 增加 量 , ne Zy, Zy = {1,2,---}, O<1 <1, 1 为 常数 , 可 以 实现 对 系统 
投放 天 敌 的 时 刻 的 控制 。Azi = r(t) rlt), Ayi = yi(t*)—-yi(t), t > 0, nT 和 (n+! 一 1)T 分 
别 表示 对 系统 采用 化 学 防治 (A Ba) 与 生物 防治 (投放 天 敌 ) 的 时 间 。 


2 ”预备 知识 与 主要 结果 


考虑 系统 (1) 的 子 系统 
Hilt) =—dizoyi(t), tA(n+l—-1)T, t#nT, 
Ayi(t) = -G+2yi(t), t=(n+l—-1)T, (2) 
Ayi(t) = pi, t=nT. 


系统 (2) 是 一 个 周期 受 迫 的 线性 系统 ， 易 得 
Di ea dit2(t-(n-DP) 


1— (1 — qi+2) se?’ 





(n-1)T<t<(n+l-1)T, 


pi(l — qi+2) se SO) 


1— (1 — qi42) - 74427 (n+1-1)T <t<nT, 





> 
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是 系统 (2) 的 周期 解 
yf (0+) = yf (nT*) = pi 


1 = (1 Gye) ee 





—dj42lT 
r= qi+2) -€ tz 


-(1- qi+2) ° 





yi ((n+1—1)*) = 区 CD = P see 
系统 (2) 以 yi(0+) > 0 为 初 值 的 解 为 
(1— gi42)"* (yi(0+) = oe + yi (t), 
(n-1)T<t<(n+l-1)T, 
yi(t) = 





Pi 一 di * 
(1— gir2)" (ye(01) Lat SC meat + y(t), 


(n+l—-1)T <t<nT. 


定理 2.1 WX (t) = (x(t), z2(t), w(t), yo(t)) 是 系统 (1) 的 解 ， 若 
d4(1 — (1 — qa)e~“47) 

















hee a2(1 + 02di + cod?) (1 — qae— dT — (1 — qa)e dT) 
ap; (1 + 01d) + c,d?)(1 — gae “4!7 — (1 ~ q3)e~*7) 1 a 
(aT ET == a = qa)e-4sT) In 1 二 | = P2maxl; 
且 
如 d4(1 — (1 — qa)e™ 7) 
a bo(1 + pode + 12d2)(1 ~ qae~44!F — (1 — qa)e~%7) 
bipi(1 + pidz + J1d2)(1 ~ gse—ds!T — (1 — gg)e7 #7 1 
(aT ipi(1 + pide + p12) (1 ~ gs as (1 — qs) ) on ) IRA 
d3(1 — (1 — q3)e 3 ) 1 一 02 
即 


P2 > max{p2maxi, P2max2} = p2max， 
则 食 饵 灭绝 周期 解 (0, 0, yilt), yo (t)) 是 全 局 渐 近 稳定 的 。 
证 明 首先 利用 系统 (1) 的 变 分 方程 来 证 明 (0,0, 寻 人 的 ,好 (的 ) 的 局 部 渐 近 稳定 性 ， 作 变 
Herit) = u(t), y(t) = u(t) + yf). FRAMAEABARA 


ur (t) ui (0) 
vat) | og | @ |, rT, 
V1 (t) vı (0) 
v2(t) v2(0) 
其 中 O(t) 满足 
dı — a1yi(t) — a243 (t) 0 0 0 
d&(t) 0 da — biyi (t) ~ bayz (t) 0 
Ae ot), (3 
dt yayi (t) ybiyi (t) —d3 j 


nazy3 (t) nbzyz (t) 0 -dz 
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县 更 (0) = TJ， 由 Floquent 理论 可 知 ， 如 果 和 矩阵 
l-u 0 0 0 100 0 
0 1- 0 0 010 0 
©= 1 E(t) 
0 0 1-g 0 0010 
0 0 0 l-4 0001 


的 所 有 特征 根 的 模 都 小 于 1， 则 周期 解 (0,0, y7 (t), yd (t)) 是 全 局 渐 近 稳定 的 。@ 的 四 个 Floquet 
HFA 


Uy = (1 ae qı) . efo (day (t)—azy3 (t))at, us = (1 D q2) ; efo (da—b1y (t) bag (t))at 
ug = (1—qs)-eS? <1, ug = (1 ga) e^ <], 
FE yt (t) 代入 积分 可 知 ， 若 满足 


da(1 — (1 ~ ga)e 7”) 
az(1 — qge~44!T — (1 — ga)e-dsT) 





p2 > 





aipi (l= gae% = (1 ~ ga)e 7) ic ) 


(ar d3(1 — (1 — q3)e-&T) i = 


da(1 — (1 ~ ga)e 7) 
bo(1 一 quae dlT — (1 e gaje744T) 





p2 > 





bipi(1 ~ Gen? = (kage) awe. ) 


(aT d3(1 — (1 = q3)e—#T) 1 一 q2 


WRTA <l. 
仿照 文献 [22] 中 定理 1 的 证 明 过 程 ， 可 以 证 明 当 pe > pomax 时 ， y(t) > ytt), y > 
y3(t)(t = 00). a 
定理 2.2 车 


a< 5, B< È A; >0, B;>0, 


A 
Az(1 — (1 ~ qa)e727) 

P2 < A A 
ag(1 — qge~42!P — (1 — ga)e™ 427) 


(a r- APU- ge AT- (1-g)e AT) a 
. Ax(1 一 (1 一 q3)e7“1T) l—-q 








人 
) = P2min1, 


B2(1 — (1 — qa)e 7") 
bo(1 — qae ` PIT — (1 — qa)e-B27) 





P2 < 


bıpı(1 ~ gae ?1!7 — (1 — q3)e ?17) 1 
By(1— (1 — q3)e-#:T) Tap 





(ar z 
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Rp P2 < min{p2min1, P2min2} = P2min? 则 系统 (1) 是 持续 生存 的 ， 其 中 


by do nb2d2 
ae, A2=d4— om 
1+ pıd2 + yids 1 + p2d2 + pods 





A, = dz 一 
yardı B =d — nazdı 
1+ 01d; + c,d?’ ae 1+ od; + cod?" 
证 明 EXE) = (xi (t), rt) yr (t), yo(t)) BAB (1) HEX (OT) > 0 的 解 ， 于 是 对 于 充分 
大 的 t， 存 在 一 个 常数 M > 0， 使 得 zi(t) < M, y(t) < M, ROR > ON, r(t) < 
M, yi(t) < M, B. 


By = dz — 








d d 
M> max{ ah 
由 系统 (1) 可 知 
ji1(t) < z1(t) (di — 21 (8), 
考虑 比较 方程 


w(t) = (di 一 Ma 的 ) t#(n+l-1)T, 

Aui(t) = -quu (t), t=(n+l—1)T, 
于 是 Zi < u(t) Hult) 一 d1， 当 t 一 00， 即 对 于 任意 的 e > 0， 当 t 充 分 大 时 zx1(t) < 
di +e, MÒR r(t) < di 十 e,t > 0; FÆ r(t) < ds+e,t>0。 令 


pi(l — git2) edt2T 














A Tague re e>0, 
考虑 以 下 方程 
y(t) > —dizays(t), t#(n+l-1)T, t¥nT, 
Ayi(t) = —q+eyi(t), t=(n+l-1)T, 
Ayi(t) = pi, t=nT. 
由 脉冲 比较 定理 可 得 yi(t) > vi(t) > yt) —e, A 4 tEAM y(t) > mi S 
yam yb (dg +€) 
Ay = d3 一 一 
Tn 
nazmı nbo(d2 + €) 
Ag = d 
2 EPE 62m, +e2M? 1+ p2(d2 +€) + p2(d2 + €)? ae 
é ybime Yoildi +€) 、0 
l+pim2+uim3 14+61(di te) +e1(d) +)? , 
bom d 
E LZ ip naz(dı +€) Ey 
1+ pome+ poems 1T 十 bo(di 十 E) 十 c2 
yamı Yb1m2 
d = 
3 TO tam? 1+pm+ pum? 20; 
di nazmı nbome 





7 1+ 69m, + com? 1+ pam2 + poms 


下 面 证 明 可 以 找到 元 ; > 0， 使 得 对 于 充分 大 的 +， 成 立 mi( > 元 ;， 分 两 步 来 证 明 。 
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步 1 下 面 将 证 明 存 在 t; > 0， 使 得 tilt) > mz， 否则 将 有 以 下 三 种 情况 : 
(i) 存在 ta > 0， 使 得 x2(t2) > mo, M x(t) < my, t>0; 

(ii) Et > 0, r(t) >m, if xo(t) < me, t>0; 

(iii) 2,(t) < mi, t > 0, 

首先 对 情况 (i) 进行 证 明 ， 选 取 充 分 小 的 cl > 0， 使 得 


as a)exp( (dı — mı — a(dz +£) — aye, 一 a2@1)T 


2 ; 
> Pidi(1 — git2.e NT _ (1 — qipa) | si 
izi Ail = (1 = qi+2)  e7XT) ? 


由 系统 (1) 的 第 三 、 四 个 方程 及 z(t) < m, T2(t) < dz +€, t>0, TË 


Hilt) < —Axy:(E), tA#(n+l-1T, tÆnT, 
Ayi(t) = -greyi(t), t= (n+l— 1)T, 
Ayi(t) = pi, t=nT. 


由 脉冲 比较 定理 得 y(t) < vilt), uM) > uF), 4t oo, HH, (t) 是 方程 


u(t) = — Aju; (t), t 天 (n +l— DT, t Æ nT, 





Av;(t) = —G+20;(t), t= (n +l 一 DT, (4) 
Avi(t) = Pi, t=nT 
的 解 ， 且 
pie™ Ni(t—(n~1)T) rT 之 I- DT 
120ga) NT? (n-1)T <t<(n+I-1)T, 
y(t) = 
i(1—qi42)-e 7 Ni(t—(n— 
ese ay (n +i- 1)T <t <n. 


AI, FET: > 0, “4t> Te H, yilt) < u(t) < v(t) +e TË 
tı (t) > zı(t)(dı ~ mı ~ a(dz +€) — a (o3 (t) ter) ~ a(v3(t) + en)), t# (n+l 1)7, (5) 
Axi (t) = =q Tı] (t), i= (n +l 一 DT: 
SN, E Z4; H (Ni+l- 1T >T, ïn >N, (Ni +1~1)T > max{T;, To} i, 对 (5) 式 
E ((n +1- 1)T, (n+ DT) 上 积分 可 得 
zı ((n +1)T) 


(n+1)T 
(dı ~ my ~ a(dz + €) — an (vt (t) + €1) ~ a2(v3(t) + ei))at) 


V 


> (ti DT) exp ( f 


(n+i—1)T 


il 


vı (n+l~ 1)T)o1, 


ai ((Ni +k +1)T) > zı ((M, + DT)ot 一 ov, 
Hk = 00, KE 3， 这 与 zi 人 的 有 界 性 矛盾 。 同 理 可 以 对 情况 (ii) 进行 证 明 ， 下 面 证 明 情 
况 (iii) 。 


第 6 期 RLH: H Monod-Haldane 功能 反应 的 四 种 群 食 饵 一 捕食 脉冲 系统 的 动力 学 分 析 1057 





$2 若 对 于 t > 五 ， 均 有 zi(b > m1， 则 取 元 7 = mi， 结论 成 立 。 否 则 ， 令 
t= inf (21t) < my}. 
1) 若是 脉冲 点 ， 即 存在 ml1 € Lys 1879 t* = (ny +l-T, QW x(t) > mi, te [t1, t*), 
且 由 x1(t) 的 连续 性 可 得 Ly (t*) 一 ?101。 选取 n2, ng E OH， 使 得 


In(e1/(M + pi)) In(e1/(M + po)) } 


一 工 
(n2 — DT > max { 一 da + yaımı + ybı -da + nagm, + nba M 


(1 一 q)opsen26T > (1 = qyapelratior > 1, 


其 中 
ô = dy — mı — a(dg +€) —a1M — azM <0. 
&T = (nz 十 nsa)T， 则 存在 ts € (t*,¢* +T], 使 得 z1(t3) > m1， 否则 考虑 (4) 式 满 足 初 始 
条 件 w(na7+) = yi(m Tt) HR, 4 


te((n-VT nT], née [ni +1,nı +n2+nz +1] 


时 ， 有 
(1 = aga)" (v(m TH) — agree) EMD +07 (0, 
(n-1)T <t<(n+1-1)7, 
vilt) = 
(1 一 Giga)? t-)) (vi(mT*) 7i noes) . @—Ai(t—n1T) 十 ux (t), 
(n+1—1)T <t<nT. 
于 是 


uilt) — vf (t)| < (M + pi) eo Nm < e, 
从 而 当 t € [(ni 十 nz 一 1)T,t* 十 TT] 时 ，(5) 式 成 立 ， 同 步 台 1 可 得 
21( + (n2 + na)T) > xi(t* + neT) 07°, 
由 系统 (1) 得 


£1 (t) > xi(t)(dy 一 ?7021 一 a(d2 +€) 一 al — aM) = ô- x1(t), tÆ (n+l-— 1)T, 
Axı (t) = —qı xı (t), t=(n+l-— 1)T. 


(6) 


对 (6) 式 在 [t*, t* + noT] 上 积分 可 得 
z1(t* + mT) > m (1 — q1)? -eT > m, 


矛盾 。 令 


t= inf {21 (¢) 2 mı}, 
WHF te ip Ent) <m, Ba(f)=m, “teei, A 


x(t) > z(t _ gq1) ts 2 eo(t-t") > mill q1)"2t™s _ es ln2+nat+l)T = Ti. 
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由 xi 人 > mi 知 ， 对 于 t > tf， 上 述 过 程 可 以 重复 。 所 以 x1(t) >m, t> te 

2) Ae RRA, Wate ftt) A, s(t) > m, Bai(t*)=m. $t e [m +l- 
DT, (+ OT), T E Ze, Htet (HOT, m(t) UL FR: 

a) Ht e (t*,(m + DT), zi(t) < mi， 类 似 于 步 2 中 1)， 可 以 证 明 存 在 ty € [a + 
DT,( 卉 十 T+ 全 ,使 得 zi1(t4) > m, > 


t= at > mi}, 


则 当 t€ ftt, e(t) <m, Hat) =m. Hte (t,t) 时 
xy (t) > m (1 — g)" -etat tl) = my, 


E r(t) >m 知 ， 对 于 上 > t 上述 过 程 可 以 重复 。 
b) 存在 ie (t*, (N17 十 中 T]， 使 得 zj1(t) > m, + 


t= inf {71(t) > mi}, 


则 当 t € (OM, a(t) <m, Half) = me Yt e (AI, (RRE, WZ 
在 (t*,t) LRT 
Zi1(t) > r(t") eC > miesnztna)T > ay, 


HÂ >m, MU4t> tf 时 ， 上 述 过 程 可 以 重复 ， 因 此 ， 当 t+ > ttf, a(t) > my. ABA 
WE, St >to, xo(t) > mz 
推论 2.1 X(t) = (z(t) r(t) y(t), yolt)) RAB) ERM, E< P, Bi > 0， 
县 pomaxl < P2 < Pomina> VW) re(t), y(t) 是 持续 生存 的 (= 1,2), Ht — œ, z(t) 一 0。 
推论 2.2 BX) = (mth rt) (4), yo(t)) BAS (1) 的 任意 解 ， 若 a < p, A > 0， 
F pomax2 < pz < pzmint， 则 Zz1(t), yi(t) 是 持续 生存 的 (i=1,2)， 且 当 # 一 00, x(t) 一 0. 


3 ”数值 模拟 与 分 析 讨 论 


对 于 系统 (1) 的 前 三 个 方程 所 构成 的 连续 系统 ， 我 们 不 能 求 出 显 式 解 ， 且 系统 (1) 不 能 化 为 
与 之 等 价 的 差分 系统 ， 从 分 析 的 角度 难以 处 理 ， 所 以 使 用 数值 模拟 的 方法 来 描述 。 考 虑 种 群 密 
度 的 非 负 性 和 和 天敌 投放 量 对 系统 (1) 动力 学 性 质 的 影响 ， 在 第 一 象限 内 作出 系统 (1) 的 分 支 图 、 
相 图 及 时 间 序 列 图 。 

取 参 数 d = 1.4, d =0.9, ds = 0.6, d4 = 0.75, a = 0.2, 8 = 0.12, al = 0.4, ag = 0.3, cl = 
0.4, co = 0.5, bi = 0.4, bz = 0.25, py = 0.5, uz = 0.4, 6, = 0.2, 02 = 0.2, pı = 0.15, p2 = 
0.3, y = 0.3, 7 = 0.5, qı = 0.02, q2 = 0.02, q3 = 0.06, q4 = 0.01, pı = 0.2, l? = 0.01, 7 =4 和 和 
PMA (1,1,1,1), WW 


Pomaxt = 13.68452135, pomax2 = 10.32658412, paminl = 10.770778565, 


Poming = 7.320205157, pamax = 13.74042805, pamin = 7.320205157. 


Wp. = 13.8， 其 相 图 与 时 间 序 列 图 如 图 1， 其 中 ， 图 1(a), (aa) 从 相 图 上 说 明了 灭绝 周期 解 的 全 
局 渐 近 稳定 性 ， 支 持 了 定理 2.1 的 结论 。 取 po = 10.5， 可 看 出 ，z2(t) 一 0， 当 上 一 oo， 其 相 图 
与 时 间 序 列 图 如 图 2， 其 中 ， 图 2(a), (aa) 从 相 图 上 说 明了 zi 人 ,的 是 持续 生存 的 且 zaz 人 的 灭 
绝 。 与 推理 2.2 吻 合 。 
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Al: 系统 (1) 的 相 图 与 时 间 序 列 图 ， 系 统 (1) 的 食 饰 灭绝 周期 解 (0,0,yi(t), 奶 (t)) 是 全 局 渐 近 稳 
Zt, (a) 相 图 (xi(t), z2(t), y(t); (aa) WB (a(t) (4), yt); (b) RF art) 的 时 间 序 列 图 ; (c) X 
于 za 人 bb 的 时 间 序 列 图 ; (d) XF y(t) 的 时 间 序 列 图 ; (e) 关于 yz(t) 的 时 间 序 列 图 























































































































































































































图 2: 系统 (1) 的 相 图 与 时 间 序 列 图 ，Z1(t), yi(t) 是 持续 生存 的 且 x2(t) 一 0， 当 一 oo (a) 相 
图 (a1 (t),x2(t), y (t)); (aa) 相 图 (x2(t), ya (t), ya(t)); (b) RF z(t) 的 时 间 序 列 图 ; (c) XF rlt) 的 时 
WESE, (d) 关 玫 (的 时 间 序 列 图 ; (e) AF yo(t) 的 时 间 序 列 图 


取 参 数 di = 5, d2 = 4.2, ds = 0.6, d4 = 0.9, a = 0.4, B = 0.28, a, = 0.71, ag = 0.51, ci = 
0.2, C2 = 0.3, bı = 0.41, ba = 0.61, AI = 0.3, H2 = 0.4, ĝi = 0.11, 82 = 0.12, pı >= 0.11, p2 = 
0.12, y = 0.95, 7 = 0.98, qı = 0.2, q2 = 0.18, gs = 0.01, g4 = 0.02, pı = 0.08, 1 = 0.1, T = 
10 和 初 值 (1, 1,1,1)， 系 统 (1) 的 分 支 图 如 图 3。 当 0 < pz < 0.4 时 ， 系 统 出 现 拟 周期 现象 ， 如 
图 4(a); 当 0.4 < pa < 6.3， 随 着 投放 量 po 的 增加 ， 系 统 (1) 出 现 典 型 的 倍 周期 分 支 到 混沌 现 
象 ， 如 图 4(b)-(e); 当 13.6 < pz < 18， 随 着 投放 量 pa 的 增加 ， 系 统 (1) 出 现 典型 的 半 周 期 分 支 
现象 ， 如 图 4(f)-(i)。 
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图 3: 
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图 4: 系统 (1) 的 相 图 ，(a) pa = 0.1， 出 现 拟 周期 现象 ; (b) po = 4, HRT AIM; (c) pa =6.1, k 
现 27 周期 解 ，(d) pz = 6.25， 出 现 4T 周期 解 ; (e) pa = 6.35， 出 现 混沌 现象 ; (f) pz = 13.1， 出 现 混沌 
现象 ; (g) po = 14， 出 现 47 周期 解 ; (h) pa = 15， 出 现 27 周期 解 ; (i) pa = 17， 出 现 全 周期 解 
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HOA AS, “Rp = 7.5 时 ， 系 统 (1) 出 现 典型 的 混沌 现象 ， 其 相 图 与 时 间 序 列 图 如 
图 5。 

目前 混沌 还 没有 统一 的 定义 ， 混 沌 的 本 质 特 征 是 对 初始 条 件 的 敏感 依赖 性 。 记 系统 (1) 过 初 
Ë (1,1,1, 1) 的 解 为 (z1(t), vo(t), y(t), yo(t))» XE (1.0001, 0.9999, 1.0001, 1.0001) 的 解 为 


( (t), X2(t), Yı (t£), Yo(t)). 
比较 两 状态 的 差异 ， 令 
e(t) = max {|x1(t) — X1(t)|, |z2(t) — Xa(t)|, l1 (6) — Y1 (0), lya(t) — OIJ, 


考虑 当 pa = 7.5 时 系统 (1) 状态 变量 对 初始 时 刻 的 敏感 性 。 由 图 6 可 知 经 过 一 段 时 间 ， 两 状 
态 变量 时 而 靠 的 很 近 ， 时 而 相差 很 大 ， 即 状态 变量 对 初始 时 刻 具 有 敏感 依赖 性 ， 也 就 是 说 
4p = 7.5 时， 系统 (1) 具有 混沌 性 质 。 
数值 模拟 佐证 了 本 文 的 结论 还 显示 出 系统 (1) 具有 复杂 的 动力 学 性 质 ， 如 随 着 参数 的 增加 ， 
系统 出 现 拟 周期 分 支 、 倍 周期 分 支 、 混 沌 和 半 周 期 分 支 等 现象 。 
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图 5: 系统 (1) 的 相 图 与 时 间 序 列 图 ，pa = 7.5, (a), (aa) 相 图 ; (b) 关于 zi(t) 的 时 间 序 列 图 ; (c) 关 
于 zz( 昌 的 时 间 序 列 图 ; (d) 关于 y(t) 的 时 间 序 列 图 ; (e) AT y(t) 的 时 间 序 列 图 


7 4 Wi 由 














图 6: Rp = 7.5， 系 统 (1) 解 对 初 值 敏感 性 ，(t elt) 
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Dynamic Analysis of a Four-Species Prey-predator Impulsive System 
with a Monod-Haldane Functional 


HUANG Wen-tao', WU Xing-jie?, LI Wei? 


(1- School of Mathematics and Computational Science, Guilin University of Electronic 
Technology, Guilin 541004; 2- Department of Mathematics, Hefei Normal University, Hefei 230061) 


Abstract: Basing on the strategy of integrated biology management and using the impulsive com- 
parison theorem, the Floquent theory and the small amplitude perturbation skill, we investigate the 
dynamical properties of a four-species prey-predator impulsive system with a Monod-Haldane func- 
tional response, impulsive ratio harvest and impulsive release. Sufficient conditions for the system of 
the two preys to be extinct and permanence are proved. Moreover, the two sufficient conditions of 
which the extinction for one of the four preys and the permanence of remaining species are given. 
Numerical simulation shows that, with the increasing of immigration, the system has more complex 
dynamics including the quasi-periodic bifurcation, periodic doubling bifurcation, chaos, periodic halv- 
ing bifurcation, etc. 

Keywords: prey-predator; impulsive comparison theorem; bifurcation; chaos; Floquent theory 
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